PURE STATES, POSITIVE MATRIX POLYNOMIALS 
AND SUMS OF HERMITIAN SQUARES 



IGOR KLEP AND MARKUS SCHWEIGHOFER 



Abstract. Let M be an archimedean quadratic module of real t X t matrix 
polynomials in n variables, and let S C R n be the set of all points where 
each element of M is positive semidefinite. Our key finding is a natural 
bijection between the set of pure states of M and S x P t_1 (R). This leads 
us to conceptual proofs of positivity certificates for matrix polynomials, 
including the recent seminal result of Hoi and Scherer: If a symmetric matrix 
polynomial is positive definite on S, then it belongs to M. We also discuss 
what happens for nonsymmetric matrix polynomials or in the absence of the 
archimedean assumption, and review some of the related classical results. 
The methods employed are both algebraic and functional analytic. 



1. Introduction 

We write N := {1, 2, . . . }, Q, M. and C for the sets of natural, rational, real and 
complex numbers, respectively. The complex numbers C always come equipped 
with the complex-conjugation involution. For any matrix A over a ring with 
involution A, we denote by A* its conjugate transpose. If A is a real matrix, 
A* is simply its transpose. Let Sym^4* x * := {A € A txt \ A = A*} be the set of 
all symmetric t x t matrices. Examples of these include hermitian squares, i.e., 
elements of the form A* A for some A 6 A txt . 

Recall that a matrix A E IR ix * is called positive semidefinite if it is symmetric 
and (Av,v) = v*Av > for all vectors v £ R , A is positive definite if it is 
positive semidefinite and invertible, and is called negative semidefinite if —A is 
positive semidefinite. For matrices A and B of the same size, we write A < B 
(respectively A -< B) to express that B — A is positive semidefinite (respectively 
positive definite). Geometrically, A £ SymlR* x * is positive semidefinite if and 
only if all of its eigenvalues are nonnegative, A is positive definite if and only if 
all of its eigenvalues are positive, and A is not negative semidefinite if and only 
if one of its eigenvalues is positive. The following algebraic characterizations 
are easy to prove: 

Proposition 1. Let A 6 SymlR* x *. 

(a) A y if and only if A is a sum of hermitian squares in M* x *; 
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(b) A j< if and only if there exist Bi, Cj G R <x< such that 

i j 

The main goal of this article is to explain how this proposition extends to 
matrix polynomials, i.e., elements of the ring R[X]* X * where M.[X] is the ring 
of polynomials in n variables X = (Xi, . . . , X n ) with coefficients from R. Note 
that in W txt , every sum of hermitian squares is of course a hermitian square. 
The reason why we speak of sums of hermitian squares in Proposition Q] is that 
this is no longer true in R[X]' X< . Note however, that A E R[X]* xt is a sum of 
hermitian squares in R[X] txt if and only if there is u G N and B & M.[X] uxt 
such that A = B*B. 

Let A be a ring with involution a t—* a* (i.e., (a + b)* = a* + b* , (ab)* = b*a* 
and a** = a for a, 6 G A) and set Sym.4 := {a € A \ a = a*}. A subset 
M C SymA is called a quadratic module in A if 

1 G M, M + M CM and a*Ma C M for all a G A 

To every G C SymR[J] ixi , we associate the set 

S G := {x G R n Vg G G : >r 0} 

and the quadratic module Mq generated by G in R[X]* X *. That is, 

M G = Oy*giPi | N G N, 5i G {1} U Cft G R[X] txt | . 

In particular, M is the set of all sums of hermitian squares in R[X]* X *. 

Given a matrix polynomial / G SymK[J] lxl and S C R n , we write / ^ on 
5 if for all x G 5, /(x) ^ 0. Likewise we use / >~ 0, / ^ 0. With this notation, 
/ G Mq implies / >z on Sq. 

In the sequel, we investigate how Mq can be used to describe matrix polyno- 
mials / G SymR[X] txt with / >- 0, / h or / j< on Sq. In Section^ the case 
G = is considered; classical results on globally positive semidefinite matrix 
polynomials in one or more variables are reviewed, and then we turn to nowhere 
negative semidefiniteness. We give a sum of hermitian squares representation 
with denominators in the one variable case and prove mostly negative results 
for the case of matrix polynomials in several variables. 

Our main results are presented in Section [3] which is devoted to the case of 
compact Sq. Actually we work under the slightly stronger assumption that 
the quadratic module Mq is archimedean (which can be enforced by possibly 
enlarging G without changing Sq). Under this assumption we describe all pure 
states on SymR[X] lxl (extremal linear forms positive with respect to Mq) as 
being of the form p i— > {p(x)v,v) for some x G Sq and v G R*. From this we 
deduce certificates, in the spirit of Proposition [H for matrix polynomials being 
nowhere negative semidefinite on Sq or positive semidefinite on Sq in the spirit 
of Proposition [TJ The latter was originally proved by Hoi and Scherer [HS] with 
entirely different methods. 
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For a study of positivity of matrix polynomials in noncommuting variables, 
we refer the reader to Helton and McCullough, see e.g. |HM| . Burgdorf, Schei- 
derer and the second author [BSS| investigate pure states and their role in 
commutative algebra. 



2. Globally positive matrix polynomials 

There are various notions of positivity for matrices. Like in Proposition [TJ 
we consider positivity of the smallest and largest eigenvalue, respectively. 

2.1. Globally positive semidefinite matrix polynomials. By Gaufi' the- 
orem, every nonnegative univariate real polynomial is a sum of two squares 
of real polynomials. The extension to univariate real matrix polynomials was 
first given by Jakubovic [Ja] and is in a different form commonly known as 
the Kalman-Jakubovic-Popov lemma [AlPj . It is one of the vast number of 
matrix factorization results obtained and used in operator and control theory 
[GKSl lGLR| IRR] . We refer the reader to [AIP] for a nice algorithmic proof; see 



also Djl 



Theorem 2 (Jakubovic). For f G SymM[Z} txt , the following are equivalent: 

(i) fyo on R; 

(ii) f is a sum of two hermitian squares in R[Z] ixi . 

Note that (ii) is equivalent to / = g*g for some g G M.[Z] 2txt , or / = YliLi v i v t 
for some v t G R[Z\* (cf. also [flnaiFHB] ). 

The multivariate version of Gaufi' theorem is Artin's solution to Hilbert's 
17th problem [MallPD] : a nonnegative multivariate real polynomial is a sum 
of squares of real rational functions. A multivariate version of Jakubovic's 
theorem (and at the same time the matrix version of Artin's theorem) was 
obtained by Gondard and Ribenboim [GRj in 1974 and reproved several times, 
e.g. IBffllFSlEENj. 

Theorem 3 (Gondard & Ribenboim). For f G SymRf-X]***, the following are 
equivalent: 

(i) / h on W n ; 

(ii) p 2 f is a sum of hermitian squares in M.[X] txt for some nonzero p G R[X]. 

Proof. From (ii) it follows that / y on {x G R n \ p(x) ^ 0} and hence 
(i). Conversely, suppose that (i) holds. By diagonalization of quadratic forms 
over a field, there exists an invertible matrix g € WL(X) txt and a diagonal matrix 
d € R(X) txt such that / = g*dg. By (i), d is positive semidefinite where defined. 
By Artin's solution to Hilbert's 17th problem, we find a nonzero p G M.[X] such 
that p 2 d is a sum of (hermitian) squares in R[X]* X *. Without loss of generality, 
we can assume that p 2 also clears the denominators in g. ■ 

In the literature cited above, one can find refinements of (ii) at the expense 
of more complicated proofs, e.g. (ii') p 2 f is a sum of squares in the commuta- 
tive ring C M[A"] ixi for some nonzero p G MLY]. Also, Gondard and 
Ribenboim [GRJ prove a bound on the number of hermitian squares needed. 
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2.2. Nowhere negative semidefinite matrix polynomials. We now turn 
to symmetric nonnegative matrix polynomials which are not negative semidef- 
inite globally, i.e., whose evaluations all have at least one positive eigenvalue. 
We start by giving an analog of Proposition [DJbj) for univariate matrix polyno- 
mials. Though this is the perfect counterpart to the well known Theorem^ it 
is to the best of our knowledge an entirely new result. 

Theorem 4. For f G Sym~R[Z] txt , the following are equivalent: 

(i) / 2< on R; 

(ii) there exist pi G M[Z] ixi such that YliPifPi — ^ is a sum of hermitian 
squares. 

Proof. It is clear that (ii) =4> (i) , cf. Proposition [TJJbJ) . 

To prove the converse, suppose first that / is diagonal, say / = 

L ' /■ 

By assumption (i), ^{-/^...,-ft} = is compact. Since we are in the univariate 
case, this implies that the quadratic module Mr f lr ..,-f t \ Q M[Z] contains all 
polynomials positive on Sr^^ ^ ^j^ |PD, Theorem 6.3.8]. In particular, with 
fo := —1, there are gij G M[Z] satisfying 

-i = E(-/*)E4- w 

i=0 j 

Observe that for each i, 

t 

fi = J2E* k fE ik eM {f} , (2) 

k=l 

where Ej k are the t x t matrix units. Thus (pQ) implies 

t t 

E /Z( Eik9i ^*f( Eik9i ^ ~ 1 = E E 9*jfi9ij ~ 1 = E 9oj (3) 
i,k=l j i=l j j 

is a sum of hermitian squares. 

Now suppose / is not necessarily diagonal. By the version of the LDU de- 
composition for matrix polynomials given in Sm2. Proposition 8], there exist 
diagonal matrices D e £ M.[Z] txt , and matrices Q G M[Z]* X *, I = 1, ...,m, 
satisfying 

(a) = C|(-/)Q, 

(b) for each iGl, -/(x) ^ if and only if for all £, D e (x) y 0. 
From (b) it follows that the diagonal matrix 

L ' -D„ _ 

tive semidefinite. If Di is the diagonal matrix with entries dj t g, j = 1, . . . , t, then 
again by |PDl Theorem 6.3.8] we deduce —1 G e \j=i t ...,t,e=i, ...,m}- Like in 
© we have G M {De} C M {jDl) ... )£ , m} for all Since M {Dli ... iDm} C M { _ /} 
by (a), we conclude as in ([3]) that — 1 G M/_fi. ■ 



is nowhere nega- 
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Unlike in the univariate case or for Proposition [UJaj) where a satisfactory 
statement on the level of multivariate matrix polynomials has been given in the 
previous subsection (see Theorem [3]) , there does not seem to exist a straight- 
forward extension of Proposition [T]Jb]) to the multivariate case. 



Example 5. Consider / G SymM[X] 2x2 . We have / ^ on W 1 if and only if 
tr f(x) > or det f(x) < for all x G R n . 

(a) Let / be diagonal. Then / ^ on W 1 if and only if there exist pi G 

is easy 
a 



R[X] 2x2 such that Y^iPifPi G 1 + -M" . Indeed, the implication ( 
(cf. Proposition [T]). For the converse implication, suppose / 



2<0 



on 



Then a > on Sr c \ and therefore p 2 a = 1 + a 



cj 
rc for some 

p G R[X] and sums of squares a, r G R[X] by Krivine's Positivstellensatz 
(see, e.g. [Mai Chapter 2] or [PD} Section 4.2]). Obviously, there exist 

diagonal K G R[X] 2x2 such that £\ h*fhi 



E 



2 

jra 


l + <7 






TC 





1 + <T 



Now 



G 1 + M, 



0- 



(b) If/ 



6 

c 

/ 



then 



+ 



/ 



tr/ 
(tr/)(det/) 



(4) 



Unable to settle the general case, we assume tr/(x) ^ for all x G R n . 
Then the diagonal matrix on the right hand side of @ is nowhere negative 
semidefinite on M n . By part (jaj) above, we obtain hi G M[X] 2x2 such that 
»://,, ( 1 • M . 



Example 6. The diagonal matrix 

~X X 
/:= X 2 

X X X 2 + 1 



G SymlR[X] 



3x3 



satisfies / ^ on M 2 and yet there do not exist pi G M[X] 2x2 with Y^iVifVi € 
1 + M . This example is inspired by [Ma, Example 7.3.2(i)] which is a modifi- 
cation of the Jacobi-Prestel example [PDj Example 6.3.1]. 

By way of contradiction, assume YliPifPi = 1 + 9 where q G M . Extract- 
ing the top left entry on both sides of this equation, we get sums of squares 
o'i,o- 2 ,a 3 ,T ^ R[^i,^2] with 

aiXi + a 2 X 2 + a 3 (X 1 X 2 + 1) = 1 + r. 

In particular, —1 lies in the quadratic module generated in Mpfi, X2] by —Xi, 
—X 2 and — X\X 2 — 1, contradicting the existence of a semiordering in R[Xi, X2] 
containing these three polynomials (cf. |Ma} Example 7.3.1]). 

We will see that a sum of hermitian squares representation with denominators 
(and weights) does exist for matrix polynomials nonnegative on a compact set 
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with archimedean corresponding quadratic module (see Subsection I3.4j) . This 
seems to mimic the situation for polynomials in noncommutative variables stud- 
ied e.g. in |KS] . where a Nirgendsnegativsemidefinitheitsstellensatz describing 
nonnegativity on a bounded set has been given, while the global case is still an 
open problem; see [KS\ Open problem 3.2] for a precise formulation. 

3. Archimedean quadratic modules of matrix polynomials 

C*-algebras A enjoy the following boundedness property: for all a E A there 
is an N E N such that N — a*a is a sum of hermitian squares (actually a 
hermitian square). In this section, we try to mimic this boundedness property 
in an algebraic context for other rings with involution. In the rigid context of 
(matrix) polynomials, sums of hermitian squares have of course to be replaced 
by a general quadratic module. 

3.1. Archimedean quadratic modules. A quadratic module M of a ring 
with involution A is said to be archimedean if 

Va£A3N£-N: N-a*aeM. (5) 

To a quadratic module M C Sym A we associate its ring of bounded elements 

H M {A) := {a E A | 3N E N : N — a*a E M}. 

A quadratic module M C Sym A is thus archimedean if and only if Hm {A) = A. 
The name ring of bounded elements is justified by the following proposition 
originally due to Vidav [Vij ; see also [Ci] for a more accessible reference: 

Proposition 7 (Vidav). Let A be a ring with involution, \ E A and M C 
Sym A a quadratic module. Then Hm(A) is a subring of A and is closed under 
the involution. 

In case A is an M-algebra, it suffices to check the archimedean condition (0) 
on a set of algebra generators. 

Lemma 8. A quadratic module M C SymM[X]' x ' is archimedean if and only 
if there exists N E N with N - £V X? E M . 

Proof. The "only if" direction is obvious. For the converse, observe that 1R[X]' X ' 
is generated as an M-algebra by X and the matrix units Eij, i,j = 1, . . . , t. By 
assumption, 

N-Xf = (N -J2 X f) + G M ' 

i i^j 

so Xj E H M 0&[X] txt ) for every j. On the other hand, E*jEij = Ejj and thus 

1 - E^j = ^2 E kk E kk € M. 

Hence by Proposition [TJ H M (^[X} txt ) = R[X] txt so M is archimedean. ■ 
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3.2. Pure states. In functional analysis, the concept of pure states is well- 
established, see e.g. }Ar} Sections 1.6 and 1.7] for the classical application to 
C*-algebras and their representations. Here we adopt these ideas to matrix 
polynomials. 

Let G C SymR[X]* x *. A linear form L : SymR[X] tx * -» R is called a state 
on (Syml[l] ixt ,M G ) if L{M G ) C R> and L(l) = 1. A state L is called pure 
if it is an extreme point of the convex set of all states, i.e., it is not a proper 
convex combination of two states other than L. 

We now come to the central result of this article. It is a matrix polynomial 
version of the well-known theorem stating that for every pure state L on a C*- 
algebra A there exists a unit vector v in a Hilbert space 7i and an irreducible 
^representation it : A — > B{TL) such that L(a) = (ir(a)v,v) for all a G A (see 
e.g. \At\ Theorem 1.6.6]). 

Theorem 9. Suppose G C SymR[X] txt and Mq is archimedean. For each pure 
state L on (SymR[X]* x *, Mq), there exists x G Sq and a unit vector v G R 
such that 

L{p) = {p(x)v, v) for all p G SymR[X]* x *. 
Proof. We extend L to C[X] txt by setting 

£(p + ig) = ~(L(p + p*) + iL(q + q*)) (6) 

for p, q G R[X]* X *. This is the unique C-linear extension of L satisfying L(f*) = 
£(/)* for all / G C[X]' X *. Let 

f ? 

M G := | >>}.'/,7<, I iV G n, 9j G {1} U G, Pj G C[X]* X * 

be the quadratic module generated by G in CLY]* X *. Then L is nonnegative on 
M%. Indeed, given / = (p + iq)*g{p + iq) with p,q £ R[X] txt and g G {1} U G, 
we have / = + q*gq) + i(p*gq — q*gp)- Applying the definition (|6|) of L, 

we obtain L(f) = L(p*gp + q*gq) G L(Mc) Q R>o, as desired. For later use 
let us observe that Mq is archimedean: write / G C[X} txt as / = f\ + i/2 
with fj G R[X]* X *. Then fj G H{Mq) C H(Mq) and i G H(M§). Hence 
Proposition H implies £T(Af§) = C[X]* X *. 

By the Cauchy-Schwarz inequality for semi-scalar products, 

J := {/ G C[X]* X * I L(/*/) = 0} (7) 

is a linear subspace of CLY]* X *. Similarly, we see that 

(p,q) := L(q*p) (8) 

defines a scalar product on C[X] txt / J , where p := p + J denotes the residue 
class of p G C[X]* X * modulo J. Let TC denote the completion of C[X] txt / J with 
respect to this scalar product. Note that Ti / {0} since 1 G" J. 

We proceed to show J is a left ideal of C[X] txt . Let / G C[X] txt . Since 
Mq is archimedean, there is some JVgN with N — f*f G Mq. Hence for all 
p G C[X] txt , we have 

0<L(p*(N-rf)p)<NL(p*p). (9) 
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This shows that L(p*f*fp) = for all p 6 J, i.e., /p G J. 
Because J is a left ideal, the map 

tt : C[X]' X * -» f^{p^Jp) (10) 

is well-defined. Here p ^ fp stands for the unique bounded linear extension 
to Ji of the left multiplication with / on C[X] txt / J , which is well-defined by 
([5]). Using the definition ([8]) of the scalar product, it is easy to see that tt is a 
homomorphism of rings with involution, i.e., a ^representation of C[X] txt on 
H. Setting w := 1 E 7Y, we observe that 

L(/) = M/K«;) (11) 

for all / G C[X] txt . 

We claim that the commutant 7r(C[X] txt )' of the image of tt in B{TL) is C. To 
see this, we take an arbitrary operator T G Tr(C[X] txt )' . Since the commutant is 
closed under the involution and T = T + r + i T ^[ , we are reduced to the case 
T = T*. By the spectral theorem, T decomposes into projections belonging 
to {T}" C 7r(C[X]* x *)'. So we can even assume T is a projection. By way of 
contradiction, assume T / and T / 1. Since T G 7r(C[X] <x< )' and w is a 
cyclic vector for tt by construction, it follows that Tw ^ and (1 — T)w ^ 0. 
This allows us to define states Li on (SymR[X] txt , M G ) by 

(ir(f)Tw,Tw) W)(l-T)w,(l-T)w) 

Ll(/) = pqp and L2(/) = iki -thi 2 

for all / G SymlR[X]* x *. One checks that L is a convex combination of L\ and 



L2. The state L being pure, we obtain L = Lj. By (jllf) . this implies 

(n(f)w,Xw) = \{n(f)w,w) = (n(f)Tw,Tw) = (T7c(f)w,Tw) = (ir(f)w,Tw) 

for all / G C[X]* X *, where A := ||Tt(;|| 2 . In particular, Tw = Xw since w 
is a cyclic vector for tt. This implies A G {0, 1} since T is a projection, a 
contradiction. 

By jLaj. Theorem 3.1], ker7r = I txt for an ideal / of C[X]. Since ker7r is 
closed under the involution, / is closed under complex conjugation. Moreover 
C[X]/I is contained in the center of (C[X]/I) txt = C[X]* x Vker tt ^ vr(C[X]* x *) 
which is C by the above. Hence C[X]/I = C and ir(C[X] txt ) ^ C txt as a C*- 
algebra. In particular, there exists x G C n such that I = {p G C[X] \ p(x) = 0}. 
Actually x G W 1 since I = I*. Also, H = ir(C[X} txt )w is finite-dimensional. 

Next we claim that tt is an irreducible ^representation. Indeed, suppose 
U is a linear subspace of Ti invariant under every 7r(/) for / G C[A]* X *. Let 
P : Ji — > U denote the orthogonal projection. It suffices to show that P G 
ir(C[X} txt y = C, i.e., tt(/)P = Ptt(/) for each / G C[X] tx *. By the standard 
trick, we reduce to the case f = f* ■ But then 

tt(/)P = Pvr(/)P = (PTr(f)Py = (tt(/)P)* = Ptt(/). 

As indicated in the commutative diagram below, tt now induces an irreducible 
*-representation 7f of C txt . This representation is unitarily equivalent to the 
identity representation 1 |Arl Corollary 2 to Theorem 1.4.4], i.e., there is a 
unitary map <3? : TC — > C* making the diagram below commute. 
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C[X] 



txt 



C[X} txt /kevn = C[X] txt /I txt 



B(H) 



n 



txt 




T h-> $T<I>* 

1 = VJ\ — Hi- &(w) 



^txt 



Let u := $(w) G C*. For each p G C[X] txt , we have 

L(p) = (tt(p)w,w) = ($7r(p)u;, &w) = (&ir(p)Q*u,u) 
= {Qit(p(x))Q*u,u} = (i(p(x))u,u) = (p(x)u,u). 



(12) 



In particular, we get tv(Auu*) = {Au,u) = L(A) G R for all A G R txt . This 
implies uu* G R* xt , that is, uu* is a real positive semidefinite rank one matrix 
and can thus be factorized as uu* = vv* for some v G R*. We can now rewrite 
([HI) as 

L(p) = (p(x)u,u) = ti(p(x)uu*) = tr(p(x)vv*) = (p(x)v,v) 

for all p G R[A~]* X *. Also note that (v,v) = L(l) = 1, i.e., v is a unit vector. 

It remains to show that x G Sg- To show this, let g G G and z£l J , Choose 
A G R* x * with z = Av. Then 



(g(x)z,z) 
since G M G . 



(g(x)Av,Av) = {A*g(x)Av,v} = L(A*gA) > 



Proposition 10. Let G C SymR[X] ixi . For eac/i x G and eac/t unit vector 
v G R*, the state L on (SymR[A > ] txt , Mq) defined by L(p) = (p(x)v,v) is pure. 

Proof. For convenience of notation, set Lq := L. Suppose there are states L\ 
and L 2 on (SymR[X]* xt , Mq) such that 2L = L\ + 1-2- Extending each Lj to 
C[A]* X * as in ([6]), this equality still holds. Define linear subspaces J, C C[X] 
as in © by J { := {/ G C[X] txt | U{f* f) = 0}. Obviously 

h n J 2 = Jo = {/ e C[X]' X * | = o}. 



(13) 

In particular, C[X] txt /J ^ C* as vector spaces. Hence Hi := C[X] txt /Ji is of 
dimension at most t. The GNS construction for Lj yields a scalar product (., .)j 



on Hi defined as in (jSj) and a ^-representation 7Tj : C[X] 



txt 



B(Hi) (cf. the 



proof of the previous theorem). Again by [La|. Theorem 3.1], there are ideals 
h C C[X] such that kervr^ = J* xi and therefore C[X] txt / ker m = (C[X]/ 'h) txt . 
In particular, 

dimB(Wi) > dim(C[X]* x Vker7rj) > t 2 dim(C[X}/ h) > t 2 
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whence dimHj > t and therefore dim 7-^ = t. Now by (|13p . we have Jo = J\ = 
J2. Therefore we have three scalar products on H := Tto = Hi = H2, and we 
find positive definite matrices G\,G2 G C txt such that 

Li{q*p) = {p,q)i = v*q(x)*Gip(x)v 

where p = p(x)v denotes the residue class of p modulo J, (this is also true for 
i = with Go being the identity matrix). Since 

v*C*GiABv = (AB,C)i = U{C*AB) = {B,~A*C)i = v*C*AGiBv 

for all A,B,C G C txt , it follows that G { A = AG { for all A G C txt , i.e., G; G C. 
More precisely, Gi = Gi(v,v) = v*GiV = Lj(l) = 1. Thus L = L\ = L 2 . ■ 

For 7^ v G R j denote by [v] the linear subspace spanned by v seen as an 
element of the real projective space P* -1 (IR) of dimension t — 1. 

Corollary 11. Suppose G C SymR[X] and Mq is archimedean. We have a 
bijection between Sq xP* _1 (1R) and the set of pure states on (SymM[X]* xt , Mq) 
well-defined by associating to each (x, [v]) G Sg x P' _1 (M), v a unit vector of 
R*, the map p 1— > (p(x)v,v). 

Proof. By Proposition [TU1 and (p(x)v,v) = {p(x)(—v),—v), the map is well- 
defined. It is surjective by Theorem [9l To show that it is injective, let 
(x, [v]), (y, [w]) G Sg x P' _1 (R) with unit vectors w,wel' satisfy 

(p(x)v,v) = {p(y)w,w} (14) 

for all p G Syml[X] (x( . Using ((HJ) with p = Eij + Eji yields vv* = ww*. Then 
[u] = [w] since (v*w)v = vv*w = ww*w = w. Setting p = Xi in dHJ), we get 
moreover = y« whence x = y. ■ 

In general, Theorem [9] and Corollary 1111 fail badly for nonarchimedean Mq- 

Example 12. Take G = 0, t = 1 and n = 2, i.e., consider pure states on 
(R[X, y], M ) where X and Y" are two variables and M is the cone of sums of 
squares of polynomials. We endow the (algebraic) dual R[X, Y] v with the weak* 
topology and consider the closed convex cone M@ C R[X, Y] v of linear forms 
L : R[X, Y] -> R with L(M ) C R> . Note that each 0/Le becomes a 
state on (R[-X", Y],Mej) after multiplication with a positive scalar (for L(l) = 
implies L = by a Cauchy-Schwarz argument). Choose a polynomial / with 
/ > on R 2 and / ^ M , for instance the Motzkin polynomial / = X 2 Y 4 + 
X 4 Y 2 - 3X 2 Y 2 + 1 [Mi Proposition 1.2.2]. Since M is closed in R[X, Y] with 
respect to the finest locally convex topology (sec, e.g., [Ma, Proposition 4.1.2(2)] 
together with [Mai Example 4.1.5]), the Hahn-Banach separation theorem [Bal 
Theorem III. 3. 4] yields Lq G M with Lo(f) < 0. As explained above, we can 
assume that Lq is a state on (R[X, Y],M ). 

Fix a double sequence (cjj)i,ieNo of cy > satisfying J^j c^LoiX 2 ^) = 1. 
Now set 

C := j L G Ml I ]T Cii L(X 2i Y^) < 1 I . 
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Then C is weak* closed since C = ^ km) {L G M% \ Eij^^i 1 ^ 2 ^) ^ *} 
and compact (for if L E C, then \L{X 2l Y 2 ^)\ < 4: for all i,j G No, and this 
implies by a Cauchy-Schwarz argument similar a priori bounds for the values 
of L on the other polynomials). In addition, both C and Mg \ C are obviously 
convex. Hence C is a cap of M@ containing Lq (see |Ph|. page 80]). 

By the Krein-Milman theorem [Ba[ Theorem III. 4.1], there exists an extreme 
point L of C such that L(f) < 0. By Choquet theory [Phi Proposition 13.1], 
L lies on an extreme ray of After normalization, L is a pure state on 

(H.[X, Y], M&). Since L(f) < and / > on M 2 , L cannot be a point evaluation. 

3.3. Positive semidefinite matrix polynomials. Now we are ready to give 
a version of Proposition [TJJaj) for matrix polynomials in the archimedean case, 
originally due to Hoi and Scherer |HS[ Corollary 1]. Using the above classi- 
fication of pure states, the proof reduces to an easy separation argument. In 
contrast to this, the original proof of Hoi and Scherer is more involved. 

Theorem 13 (Hoi & Scherer). Suppose G U {/} C Syml[l] ixi and M G is 
archimedean. If f >- on Sq, then f G Mg. 

Proof. In the terminology of Barvinok |Bal Definition III. 1.6], Proposition \7\ 
together with the identity 4s = (s + l) 2 — (s — l) 2 shows that Mq is an archi- 
medean quadratic module if and only if 1 is an algebraic interior point of the 
convex cone Mq C SymR[X] <x< . Recall: / is an algebraic interior point of Mq 
if for every p G Syml[J] ixi there exists e > with / + ep G Mq. 

Suppose / ^ Mq. We will find x G Sq such that f(x) )f 0. The existence 
of an algebraic interior point of Mq allows us to separate the convex sets Mq 
and K>o/ by the Eidelheit-Kakutani separation theorem |Bal Theorem III. 1.7]. 
More precisely, there exists a state L on (SymMfA 5 "]^*, Mq) with L(f) < 0. The 
set of all such states is weak* compact by Tikhonov's theorem (cf. the proof of 
Alaoglu's theorem \B&\ Theorem III.2.9]). Hence by the Krein-Milman theorem 
[Ba, Theorem III.4.1], L can be chosen to be pure. 



By Theorem [9l there exists x G Sq and a unit vector v G Mr such that 
L(p) = (p(x)v,v) for all p G SymM[A > ]' x '. In particular, (f(x)v,v) = L(f) < 



Corollary 14. Suppose G C SymM[X]* x * and Mq is archimedean. For f G 
SymR[X] <x< , the following are equivalent: 

(i) / h on S G ; 

(ii) / + e G M G for all e G M> . 

For t = 1, Theorem [13] specializes to Putinar's Positivstellensatz [Puj . To 
avoid possible confusion, we use the letter Q to denote quadratic modules in 
commutative rings with trivial involution. For instance, if G C M[X] we denote 
the quadratic module generated by G in M[X] by Qq, i.e., 



as desired. 




Note that Q G = M G for t = 1 but Q G C M G for t > 1. 
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Corollary 15 (Putinar). Suppose Gu{/} C R[X] and Qg is archimedean. If 
f > on Sg, then f £ Qg- 

Clearly, if Qg is archimedean then Sg is compact. The converse is false even 
for finite G Q R[X] as shown by the Jacobi-Prestel example |PDl Example 
6.3.1]. Nevertheless, there is an intimate connection between compactness and 
the archimedean property established by Schmiidgen |Sml| . To describe his 
result, we introduce the following notation: Given a set G = {g\, . . . ,g m } f= 
R[X] of m distinct polynomials, let G = {g{ x • • • I / 5 € {0, l} m } denote 
the set of the 2 m — 1 nontrivial products of the gi . 

Theorem 16 (Schmiidgen). Suppose G C M.[X] is finite. Then Sq is compact 
if and only if the (multiplicative) quadratic module Qq is archimedean. 

As an important special case, we obtain that for a singleton G = {g} C R[-X], 
Sg is compact if and only if Qg is archimedean. This continues to hold if G 
has exactly two elements |PD| Corollary 6.3.7]. For this and other nontrivial 
strengthenings of Schmiidgen's theorem due to Jacobi and Prestel we refer to 
[PDj Chapter 6]. These results allow us to deduce that Mq is archimedean in 
such cases by the following proposition. 

Proposition 17. If G C RLY], then Qg is archimedean if and only if Mq is 
archimedean. 

Proof. To prove the nontrivial direction, suppose that Mq is archimedean, i.e., 
N ~ T,i=i X f = EjPjPjdj for some N G N, Pj 6 R[X] txt and gj G G U {1}. 
Since the trace of a hermitian square p*pj is a sum of squares in M[X], it 

follows that N — J27=i X? = \ Ylj^ T (PjPj)Sj 6 Qg- Hence Qg is archimedean 
by Lemma [HI ■ 

There does not seem to exist a viable generalization of Schmiidgen's theorem 
for general finite G C SymM[X] txt . It does not make sense to consider G 
because products of positive semidefinite matrices are not symmetric in general, 
let alone positive semidefinite. If G is a singleton, we have G = G but still Sg 
compact does not imply Mq archimedean. 

Example 18. Let / € SymR[A > ] 3x3 be the diagonal matrix from Example [6j 
Then 5{_/} = is compact but M|_j| is not archimedean. Otherwise Theorem 
1131 would imply —1 € Mrf\ which is not the case as seen in Example El 

We now briefly turn to positivity of not necessarily symmetric matrix poly- 
nomials. For this we need the following classical lemma [Br|, Section 6.3] (see 
also [Sw\ Theorem 5.3]). 

Lemma 19 (Brumfiel). Let 1Z be a commutative Q-algebra and Q C 1Z a 
quadratic module. Then Hq(JZ) is integrally closed in TZ. 

Let G C R[X] and / G R[X] <X< . The quadratic module generated by G in the 
commutative ring R[X, f] endowed with the trivial involution will be denoted 
by Q f G . Observe that Q G C M g if and only if / = /*. 
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Theorem 20. Suppose G C R[X], Qg is archimedean and f G RLYf. If for 
all x G Sg, all real eigenvalues of f(x) are positive, then f G Q G . 

Proof. Let qj G R(X)[y] be the minimal polynomial of the matrix /. Note 
that qj G M[X,Y"] by Gaufi' lemma since qj divides the (monic) characteristic 
polynomial of / by the Cayley-Hamilton theorem. 
Now 

Y > on {(x, y) G S G x R [ q f (x, y) = 0} = S GU {q f ,-q f }- 

We claim that QGu{q f ,-q f } = Qg + ^[-^> ^ s archimedean, or equivalently, 
the quadratic module Q G is archimedean in R[X,f] = R[X ,Y]/(qf). Indeed, 
since / is integral over R[X] and H qS [R[X,f]) D # Qg (M[X]) = R[X] is inte- 
grally closed, we have H nf (R[X, /]) = R[X, /]. 
By Corollary HS1 V G Q G u{ 9/ , 

Plugging in / for Y yields / G Q G . ■ 

Corollary 21. Suppose G C R[X] and Qg *s archimedean. For f G R[X]* X *, 
i/ie following are equivalent: 

(i) /or a// x G Sg, aii real eigenvalues of f(x) are nonnegative; 

(ii) f + e e Q f G for all e G R> . 

Proof. (i)=^(ii) follows from Theoreml20[ For the converse, it suffices to observe 
that for p G Q G , all real eigenvalues of p(x) are nonnegative for all x G Sg- 
Indeed, suppose p = Y^j hj(X ,p) 2 gj for finitely many gj G G U {1} and hj G 
R[X, Y]. Let A G R, + v G M* and p(a?)« = Au. Then 

Xv = p{x)v = ^2hj(x,p(x)f 'gj(x)v = ^hj{x, \fgj{x)v 

j 3 

whence A = ^ • hj(x, \) 2 gj(x) > 0. ■ 

3.4. Matrix polynomials not negative semidefinite. We conclude this ar- 
ticle with an application of Theorem [13] yielding a version of Proposition [Djb]) 
for matrix polynomials. 

Corollary 22 (Matrizenpolynomnirgendsnegativsemidefinitheitsstellensatz) . 
Suppose G C SymRLY] txt and Mq is archimedean. For f G SymR[X] txt , the 
following are equivalent: 

(i) / 2< on S G ; 

(ii) there exist pi G RLY] ixi such that 

J>i/Pi€l + M G . 

j 

Proof. (ii)=>-(i) is immediate from Proposition [Hjb]) . For the converse, note that 
/ ^ on Sg if and only if SGu{-f} = - m this case, — 1 >- on SGu{-f}- 
Since Mq and therefore Mgu{-/} is archimedean, Theorem [13] implies that 
— 1 G Mgu{-/} which is exactly what we need. ■ 

Acknowledgments. We thank Ronan Quarez for interesting discussions which 
led us to the discovery of Theorem 4. We also appreciate the anonymous referee 
for his careful reading of the manuscript and his corrections. 
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